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ABSTRACT 



This thesis addresses the problem of estimating the distance which 
must be traveled to visit all the targets in a search area. The search 
area is assumed to be a compact area drawn from an infinite Poisson 
field. Targets are visited in "nearest uninvestigated neighbor" order. 
A simplified analytical model is developed for the unbounded 
investigation. This model is used to provide heuristic justification 
for simulation results obtained for the bounded investigation. 
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I. INTRODUCTION 



A. GENERAL DESCRIPTION OF THE PROBLEM 

This thesis addressed the problem of estimating the distance 
which must be traveled to visit targets which are randomly distributed 
over a large area. The problem is interesting because many practical 
search operations require the investigation of targets after detection. 
It is often necessary to estimate the time which will be required for 
the entire operation before the targets are even detected. 

Classical search theory addresses cases where it can be assumed 
that there are few targets in the area, and that all are of equal 
interest when detected. Later classification may rule out some targets 
as false targets. This is not usually obvious to the searcher. Any 
detection will stop the search. Subsequent classification as a false 
alarm may result in its restart. 

A variation on this case occurs when there are many false targets, 
but they may be classified immediately upon detection. If the target 
may be immediately classified as false, then the search can proceed 
unimpeded. False targets can be treated as if they had never been 
detected. If all false targets are immediately classified as false, 
then this case becomes equivalent to the classical case of sparse 
targets. The search continues until the true target is found or the 
area is exhausted. When the first target is detected which cannot be 
dismissed as false, the search is stopped and that target is pursued. 

There is a third case which shall be addressed. It is inter- 
mediate between the first two. It can be defined as follows: targets 
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can be easily classified, but only at very close range. The searcher 
can detect targets at long range, but, for classification, must take 
a "close look" at which time he knows immediately whether or not he 
has found the object of his search. If he has, he stops, if he hasn't, 
he keeps searching. If the targets are dense enough, or the searcher's 
detection range great enough, the searcher may always have uninvestigated 
targets in sight. He may not really "search" at all in the sense of 
classical theory, he may instead travel directly from target to target 
without any regular "search pattern." The problem is now one of 
"investigation theory" as first defined by Nuendorfer [ll. The 
investigation will continue until the desired object has been found, 
or until all of the targets have been visited. In either event, the 
investigator's track will have gone to the location of each prior 
false alarm. 

This is the crux of the problem. If the positions of the targets 
are unknown until detection, then the length of the entire path will 
be uncertain. If the target positions are random variables, then the 
length of the path required to investigate the targets will also be 
a random variable. 

B. EXAMPLE 

The case described above is typical of many Coast Guard operations. 
The Coast Guard is often required to search for a particular boat. 

Usually the boat can be positively identified only by name or registra- 
tion number. Small boats which have been reported overdue, and boats 
suspected of law enforcement violations are two classes which must be 
searched for in this manner. There may be many boats which fit a 
general description, and the searcher may have to check out all of them. 
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A simple scenario can be constructed which illustrates the problem 
Consider a Coast Guard Cutter with embarked helicopter assigned to 
patrol a ten thousand square mile area for fishing violations. The 
helicopter must be fueled tonight for tomorrow's first light patrol. 
The aircraft commander must decide how much fuel to load. Too little 
fuel means an incomplete patrol. Too much fuel means that something 
else must be left behind. Experience has shown that rarely are more 
than twelve or fewer than eight ships seen. Ten is the average number 
The aircraft commander has no way of predicting the future position 
of any of the fishing vessels, but he must still make an important 
operational decision, based on the length of the path which he expects 
to fly. The result developed in this paper would be directly 
applicable in this scenario. 
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II. ANALYSIS 



A. ANALYTICAL FORMULATION OF THE PROBLEM 

Two basic questions must be addressed first off: 

1. What distribution is to be assumed for the positions of the 
targets? 

2. In what order are the targets to be visited? 

In this thesis it is assumed that the area considered is a portion 
of an infinite Poisson field with targets distributed with spatial 
Poisson density p. This means that the number of targets in any 
finite area. A, will be a random variable, with a Poisson distribution 
and mean pA. P.A.W. Lewis and 6.S. Shedler [2] have shown that the 
targets occurring in an area as the result of this process will be 
uniformly distributed within the area. This is to say that if the 
entire area is divided up into many small boxes, each large enough to 
hold only one target, then each box has an equal chance of actually 
holding one target. This distributional assumption is attractive 
because it assumes the minimum amount of information about the target's 
location. In this sense it is the "worst case"; more information can 
only improve the situation. In actual Coast Guard practice, there are 
many searches conducted when uniformity is the best assumption which 
can be made; the actual distribution of targets, including false targets, 
cannot be distinguished from what would have been expected from a 
uniform distribution. 

A related question is that of target motion. In this thesis it 
will be assumed that the targets move so much more slowly than the 
investigator that they may be considered stationary. 
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The question of the order in which targets are to be visited is 
more complicated. If the targets are stationary, there is one unique 
optimum order resulting in the minimum total path length. The problem 
of determining this optimum ordering is called the "travelling 
salesman problem." The travelling salesman solution, while ideal, 
is not practical for large problems. There is no efficient method of 
solution. Every current method is fundamentally equivalent to trying 
all possible orders [ 3 ]. There are some shortcuts for small numbers 
of points; the fact that points on the perimeter will be taken in the 
relative order that they are found around the perimeter is an example. 
But these shortcuts do not work for large numbers of points with 
complex geometry. Determining the optimum path through twenty points 
will require a tremendous amount of calculation. 

This thesis employs the myopic decision rule in which the investi- 
gator will always go to the nearest point which has not yet been 
visited. This decision rule is justified on two grounds: 

1. It is easily employed in actual searches; the searcher need 
only remember what targets have been investigated, and then select the 
closest of the remainder. 

2. Balut [4] has shown simulation results indicating that the 
path resulting form the myopic decision rule is only about 15% longer 
than the optimal path for investigations where the investigator is 
much faster than the targets. This procedure is often used in actual 
searches . 

Finally, the assumption is made that the investigator's detection 
range is much greater than the inter-target spacing. This means that 
the investigator will always know the location of the nearest uninvesti- 
gated target. 
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B. MODEL FOR AN UNBOUNDED MYOPIC INVESTIGATION 



Because the area considered is a portion of an infinite Poisson 
field, it is profitable to first consider the properties of an investi- 
gation conducted in an infinite field. The investigator will start 
at one of the points in the field and then visit its closest neighbor, 
then the closest neighbor remaining, and so on. The investigator will 
always travel to the nearest uninvestigated target (n.u.t.). Real 
searches never detect an infinite number of targets, but this is not 
a bad approximation to the case when the area is very large and there 
are many more targets detected than can be investigated in the course 
of a mission. 

Consider a Poisson field with spatial density p. Call the set of 
points I = { p Q , p-| , P 2 , P 3 , ..., P n > an investigation starting at p Q . If 
P n = Pg, i.e. the investigator returns to the point where he started, 
then the investigation is a tour. The investigation is conducted 
according to the following rule: the investigator will always visit 
the closest point to which he has not yet been. Let r. be the distance 

k n 



from p i _ 1 to p.. Let £ k = l 



r i 



Let t 



n 



= l 

i=l 



r i 



when 



Pn = 



z, is the cumulative distance convered in visitina the first k ooints; 
k 

t is the length of a tour visiting all the points and returning to 
its starting point. 

The expected length of the first radius, r^, may be calculated as 
follows. The probability of a target being found in a small area, dA, is 

P(dA) = pdA. 

Thus the probability that a target will not be found is 

P(dA) =1 - pdA. 



11 



Let F(A) be the probability that there is no point in some large area A. 
The probability that there is no point in (A + dA) is 

F(A) • P(dA) = F(A)(1 - pdA). 

But this is also F(A + dA) 

F(A)(1 - pdA) = F(A + dA) 

thus 



F < A + - F(A > - -PF(A). 

As dA 0, then this becomes 



dF(A) 

dA 



— pF (A) 



so that 

F(A) = ce' pA 

where c is the arbitrary constant of integration. 

As F(0) = 1, then c = 1 and 

F(A) = e' pA . 

The area considered so far has been perfectly general, of arbitrary 
shape. 

Now consider A to be the area circumscribed by a circle of radius r, 

2 

A = tt r . Let A-j be the smallest such circle tangent to an adjacent 
point. Let r-j be its radius. The tail distribution of A.| , F^ (t) is 
defined as follows: 

F A (t) = P(A ] > t) 

2 

Let t = ttx , where it is the area circumscribed by a circle of radius x. 
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Then 



Define 



But 



The expected 
Integrating 

E h 



f a (t) 
rt i 



= P(A 1 > t) 

= P(iTr 1 2 > TTX 2 ) 

= P^ 2 > x 2 ) 

= P(r-| > x) as x,r^ >_ 0 . 



G (x) = P(r, > x) 
r l 1 

= F /\ * 

A 1 



F a (t) = e pt 
A 1 2 
= e " PlTX 



G (x) = e” pTrX . 
r l 

value of r^, E(r^), is defined as 

-dG , > 

r-, (u) 

> = I [—as — ] u du 



■f 



oy parts yields 
) = | G,. (u)du 



c °° r °° 2 

= / G„ (u)du = / exp [-Pttu ] du 

JO r i Jo 



1 



2/P 



where G (x) = P(r, > x) 
r l 1 



p = Spatial Poisson density. 
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Thus the expected radius of the smallest circle containing an adjacent 
point is 



E(r 1 ) = 



1 

2/^p 







Figure I I -1 

The expected values of successive jumps are much more difficult 
to determine. The distribution in the fresh (unsearched) area is still 

F(A) = e- pA 

but the new area is no longer a circle. Let x-| be the actual distance 

between p Q and its nearest neighbor, p^ . Then a circle of radius x-| , 

centered at Pg, is known to contain no other points besides Pg. The 

new point, p^ , lies on to this circle. Area included in this circle 

"doesn't count" in the search for new points and must be ignored in 

calculating the distribution function for the next interpoint distance, 

G" (x). Let H.(x) be the function whose value is the new (uninvestigated) 
r 2 1 

area lying within a circle of radius x, centered at P. . 
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Here 

A 




Figure I I -2 



r 2 < 2r^ , thus 

2 

2 ~~ H.j (r^) — ir{r^ 



i 

- r, ) 2 - {2r, (-^-y - 1 ) cos' 1 
1 2r ] c 





15 



Thus 



2 _ -pH n (x) 

H n (x) = irx and G (x) = e u 
u r ] 

It can easily be seen that 

H-j (x) = ir(x^ - x-j 2 ) if (x >_ 2x^ ) 



? ? 2 , / 2 
= ^ - x, ) - {2x 1 (-^y)cos" 1 ^- - 757- A - - 2 -x- } 

1 1 2x^ ^ X 1 ^ X 1 4x-j d 

if (x < 2x-| ) where x^ is the actual observed distance between and p.| 

Let G (x) = P(r 9 > x). 
r 2 ^ 



Then, as before 
G r (x) = e 



•pH, (x) 



and 



E(r 2 ) = 



fo 

• 2x 



-pH, ( u) 
e du 



1 



exp(-p{ir(x - ) - [2x-| ~ 1) 



cos 



-1 x 



2x, 



2x^ “ 2x-| 



r. 



4 x n 



]}}dx 



2 2 

+ / exp { -pir(x - x, ) }dx 

J2x 1 1 

This expression is difficult to evaluate analytically, and the expected 
values of r^, r^ ... are even more difficult to obtain. This is because 
the expected value for r^ depends not only upon r-j and r 2 ... r^_-| , 
but also upon the angles between successive radii. An investigation 
which doubled back on itself, Figure II— 3 (a) , would have a different next 

expected radius than one which did not. 
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Figure 1 1 -3 

While the total sanitized area is the same in both cases, it is 
not distributed according to a common function of radius. 

Because the distributions and moments are so difficult, approxima 
tions would be useful. An argument can be advanced to support the 
conjecture that E(r^), i = 1, 2, 3, ... are asymptotically linear as 
i ->■ oo. Because new points will be located in an uninvestigated area, 
the process will tend to creep away from sanitized areas and toward 
unsanitized areas. This will usually put the investigator at the 
boundary of a fresh area, with a constantly expanding sanitized area 
behind him. Asymptotic linearity may be shown easily in a restricted 
model based on this behavior. Consider a completely deterministic 
process where points appear on a straight line. The spacing between 
the points will be determined as follows: 
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1. will be arbitrary, 

2. r. will be such that 



irr 



2 ■ - t,(r. 2 



i-1 



) - 



r. 



{2 r. ( — j - 1 )cos 

2r i-l 



•1 r i 



2 Vl 



r . 

^7 





4 r i-i 



} 



for all in i > 1 . 

This is equivalent to circumscribing a circle of radius r^ about the 

first point, and then adding other points so that the new area added 

to the process at each step is equal to the initial area enclosed in 

the first circle. The process can be visualized as a sequence of 

circles proceeding down a line. Each new circle is centered on the 

intersection of its predecessor and the centerline. Each circle is 

sized so that its area, less the area of its predecessors, is equal to 
2 

r-j . Refer to Figure 1 1 -4. 




Figure I I -4 

Numerical calculation for all initial values selected between 
-1 50 

r-j = 1 x 10 and r^ = 1.0 x 10 have shown that the sequence 

r l ’ r 2’ r 3’ r n ’ ••• 1S asymptotic to 1.2814 r-j . 
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While one would never expect to see an actual realization develop 
along a perfectly straight line, an extensive computer simulation by 
Laidlaw [5] does show that many investigations do develop in a way that 
could be characterized as random deviations on either side of a mean 
line of progress. It seems reasonable that, if this were so, and the 
investigation did proceed as variations about a mean line, then the 
actual limiting value might be r n -*■ cr.| , where c >_ 1.2814. Laidlaw' s 
simulation yielded c-j = 1.405 based on a myopic investigation conducted 
in a finite plane containing 10,000 points, with the investigation 
stopping when the boundary of the plane was closer than the nearest 
target, and with the sweep width much greater than the expected target 
spacing. 

The assumptions for this result were: 

1. The points visited are located in an infinite Poisson field. 

2. Detection range is much greater than target spacing. 

3. Myopic search is employed. 

Under these assumptions the length of the total path appears to 
be asymptotically convergent to 

l - ncr, 
n 1 

_ n(l .405) 

2/p 

The model is asymptotically linear in n and is called the asymptotic 
linear model of the unbounded investigation. 
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C. APPLICATION TO THE BOUNDED INVESTIGATION 

The linear asymptotic model cannot be directly applied to the 
bounded investigation. The linear model requires that the spatial 
Poisson paramenter, p, remain constant in the unsearched area. This 
will not be true in a finite area, because the finite area will 
contain only a finite number of points. 

Nevertheless, the linear asymptotic model does suggest an indirect 
approach. Consider a bounded investigation to consist of two parts: 

1. An investigation of the points near the center of the area 
where the unbounded linear model may be applied. 

2. An investigation of the remaining points near the boundary. 

It seems reasonable to suppose that the bounded investigation might 
first behave like the unbounded case, when there are a large number 
of points from which to draw. After a while the points will have 
been "thinned out" considerably, and the investigator's path will take 
him back to "mop up" the remainder. 

A simulation was constructed to test this conjecture. A number 
of points were generated in a square. An investigation was started 
at a point near the center. The path took the investigator to all 
the points in the area, and then back to his initial point. Such a 
closed path is called a tour. Plots of il^/kE ( r -| ) - vs - k for 
k = 1 , 2, 3 . . . n, n = 10, 30, 60 shown in Figure III-5. A cumulative 
plot of t n /nE(r-|) - vs - n is shown in Figure III -6 . 

Two things can be seen in these graphs: 

1. The "plateau" in the center of the curves in Figure I I I -5 becomes 
more pronounced as more points are added to the area. This is a 
reasonable result. 
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KEY: O+N-20, +-*tf=40, »-*^=60 



WI0 = {L„)*(K*EIR .]) 
2.04 * 



1.30 



1.56 



1.32 



1 . 08 i 



0 . 84 




o o 

+ o 



♦ o<3 



0 10 20 30 40 50 60 

K, THE HUMBER OF POINTS VISITED SO FAR. K = 1,2,...,V 



RATIO OF CUMULATIVE PATH LENGTH, L , TO (K*ElR ,]) -VS-K,K=1,2, . . . ,N 



Figure I I -5 
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I O 

i i i ; i i i i i i i i i 

0 10 20 30 40 50 60 

/V, THE TOTAL HUMBER OF POINTS IN THE PLANE 

RATIO OF TOTAL TOUR LENGTH , T^, TO (N^ElR / l)-VS^ N 



Figure 1 1 -6 
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2. The final cumulative length is approximately 

t p = 1 .95 ( 4p = 1 .95nE(r 1 ) , n > 7.' 

Plots of repeated simulation trials suggest that t p is normally 
distributed with mean 

if = 1 95 mA 
z n 2 

S f 2 = 0.06 
n 

Stemleaf plots for the cases n = 20, 40, and 60 are shown in 
Figure I I 1-7 . 

This is the expression for ¥ when n is known. Generally n will 
not be known, but the assumption that the search area is taken from 
a Poisson field will enable the calculation of E(t n ) from p and A. 
The expected number of points occurring in an area A is E(n) = pA. 
Thus 

E(tjn) = l.95uE(r 1 ) 

1 



= 1.95 n 



2/7 



and 



oo 

E(t ) = l 1 .95 k — 1— f (k) 
k =0 2/7 n 



= 1.95 -4= l k f (k) 
2 / p k =0 n 



where f n ( k ) = P(n = k) 



= 1.95 



2/7 



E(n) 



= 1 .95 — — • pA = 1 .95 

2/7 2 
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Stemleaf Plots 



nE( r 1 ) 



for n = 20, 



40, 



60 



Figure I I -7 
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Similarly, 



2 n /\^/2 

S, « 0.06 



III. EXAMPLE 



The procedure developed may be applied to the example given in 
Section II. 

Assume that the search area is a square area, 100 miles on a side, 
with a spatial Poisson parameter of 

p = 1.1 x 10' 3 . 

The expected length of a tour will be given by E(t n ) = 1.95 
thus 

E(t ) = 323 miles, 
n 

To check this result, eleven points were selected at random, with 
one being the searcher. The positions, measured from the southwest 



corner, were: 






Index 


x_ 


£ 


1 


32.1 


51.1 


2 


60.5 


00.4 


3 


90.5 


17.3 


4 


49.8 


18.2 


5 


18.4 


99.2 


6 


65.3 


72.9 


7 


40.8 


12.8 


8 


51 .6 


22.2 


9 


69.7 


99.3 


10 


58.0 


74.0 


0 


18.0 


19.8. (searcher) 
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These points are plotted in Figure III-l. It is assumed that the 
searcher's detection range is greater than the greatest inter-target 
spacing. 

The myopic investigation starting at point "0" is shown by the 
solid line in FigureIII-2. I = (0, 7, 4, 8, 2, 3, 6, 10, 9, 5, 1, 0}. 
The actual length of the myopic investigation starting at pt 0 is 
328.66. In this case the myopic path is also the optimal path through 
all the points, and 328.66 is the minimum length. In this case n = 11 
and 

S 2 * .06 Jp 1650 
t n 

S » 40.62 
r n 

So the actual value is well within one standard deviation of the 
predicted value. 
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9 



10 



'6 



0 

Searcher 



* 8 

4 



7 



.2 



Figure III-l 

Point Locations for Example Problem 
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Figure I I I -2 

Myopic Tour Through Example Problem 
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IV. SUMMARY AND CONCLUSIONS 



This thesis has presented a method of estimating the length of a 
path required to visit all of the targets in a search area under three 
conditions : 

1. The targets are distributed uniformly within the area. 

2. The "nearest uninvestigated target" (myopic) decision rule is 

used. 

3. Detection radius is much greater than inter-target spacing. 
Simulations results have shown that the expected distance of a complete 
tour of all points under these assumptions is predicted very well by 

a simple function of the area and number of targets, if that number 

is known, or by the area and average target density, if the number of 

targets is not known. 
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[5] A CENTER OF EACH POINT SET. CUMULATIVE AVERAGE RESULTS ARE OUTPUT T 

d ROUGH 

[6] A ’ CUMS . OFFLINE PRINTING OF INTERMEDIATE AND FINAL RESULTS IS COST 

ROLLED 

[7] A THROUGH THE VARIABLE OPRINT. 

[8] SIZE-NPOINTS- 1 

[9] RATIOS*- 0 

[10] CMS ' CLOSIO PRINTER OFF ' 

[11] ' 

T 

[12] * NUMBER OF TRIALS THIS RUN ' ; NT RIALS 

[13] ' NUMBER OF POINTS GENERATED PER TRIAL NPOINTS 

[14] 'SPATIAL POISSON CONSTANT, LAMBDA ' •,LAMBDA—NP0INTS*1Q*& 

[15] EX LENGTH^ 2 * LAMBDA * 0 . 5 

[15] 'EXPECTED LENGTH OF FIRST JUMP ' ; EXLENGTH 
[17] » 

T 

[13] CUMS* 0 

[19] IN DEX-NT RIALS 

[20] STEP : NP-N POINTS 

[21] B-(NPOINTS?100Q),Ll.5](NPOINTS?100Q) 

[22] COUNT-lipB 

[23] STEP2 :A—,BtCOUNT ; ] 

[24] -+STEP3x\(, 5 OS | 500-i4[l] )*50>| 500-i4[2] 

[25] -STSP2*\C0UNT+C0UNT-1 

[26] -STEP 

[27] STEP3 : 3*- 1 0 4 (C0UNT~l)e3 

[28] RATI OS-RATIOS, A NUT 2 B 

[29] CUMS+CUMS+OUTPUT , [ 2 ] ( OUTPUTl ;6]*2),[1.5]( OUTPUTl ; 7 ] *2 ) 

[ 30 ] -STEP* \INDEX— INDEX - 1 

[31] CUMS— CUMS * NT RIALS 

[32] CUMSL •,8]-(NTRIALS*NTRIALS~l)*(CUMSl ;8]-0/MS[ ;6]*2) 

[33] CUMS [ ; 9 ] ■ -{ NT RIALS iNT RIALS - 1 ) * ( CUMS [ ; 9 ] •CUMS [ ; 7 ] * 2 ) 

[34] CH-NUT2CH, ' ,L V AR,R VAR' 

[35] Q-A+TRIALSFMT FMT CUMS 

[36] -N OFF LINE * x OPRINT - 0 

[37] OPRINT A 

[38] NOFFLINE : CMS 'CLOSIO PRINTER ON' 

[39] RATIOS-1 + RATIOS 
7 
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NNUT 2[Q]V 

V PATHLTH-ST ART STOP RUT 2 POINTS ; DIST ; DIST2 •, UNVISITED ; RE XT ; ROW ; RUM 
NUM21 ;START iSTOP 

[1] A THIS FUNCTION FINDS THE MYOPIC INVESTIGATION THROUGH THE 

[2] A POINT SET ' POINTS' STARTING AT 2+ ST ART STOP AND ENDING AT 

[3] '2tS7V3rtTSTOP. OUTPUT IS THROUGH THE ARRAY ' OUTPUT' AND 

[4] A CONSISTS OF THE LIST OF POINT INDICES ( FROM 'POINTS'), POINT 

[5] A COORDINATES, JUMP LENGTHS, CUMULATIVE LENGTHS, AND RATIOS 

[6] A OF CUMULATIVE LENGTH TO (N*EXLENGTH) , WHERE 'EXLENGTH' IS THE 

[7] A EXPECTED LENGTH OF THE FIRST JUMP . 

[8] HUM-SIZE 

[9] NUT JUMPS - 0 

[10] NUT LENGTHS - 0 

[11] LENGTH-0 

[12] NUTPATH— NOW-1 

[13] LIST2-(NUM+l)*LIST-\N-2+l*pPOINTS 

[ 14 ] ST ART-2* ST ART STOP 

[15] STOP-' 2* ST ART ST OP 

[15] DIST-DI STANCE ST ART, ill POINTS, ill STOP 
[17] DIST— DIST -ZAP* ( i# ) ° . = ( \N ) 

[13] UNVISITED-0,((N’2)ol),0 

[19] STEP :DIST2-UNVISITED/ DIST 

[20] NEXT-l*i,DIST2iNON-,l 

i 21 1 LENGTH-LENGTH+DUMMY- , DIST2i NOW ; NEXT ] 

[22] NUT JUMPS-NUT JUMPS , DUMMY 

[ 23 ] NUTLENGTHS-NUT LENGTHS , LENGTH 

[ 24 ] NOW- 1 1 ( NEXT - 1 ) $U INISITED/ LIST 

[ 25 ] UNV I SIT EDiNOW 1-0 

[26] NUTPATH-NUTPATH , NOW 

[27] -STEP * \ 0 *NUM-NUM~1 

[28] NUTPOINTS-START , ill POINT Si 1 4 NUTPATH ~ 1 ; ] , [ 1 ] STOP 
i 29 ] NUTPATH-(NUTPATH-l ) , 0 

[30] NUTLENGTHS-NUT LENGTHS , LENGTH+PATHLTH- , DISTiNOW ; N ] 

[31] NUT JUMPS-NUT JUMPS , PATHLTH 

[32] OUTPUT-i LIST2 , 1 ) , [ 2 ] NUTPATH , [ 2 ] NUTP0INTS,i2l NUT JUMPS, [2] 

NUT LENGTHS , [ 1 . 5 ] ( NU T LENGTHS * EXLENGTH * ( \ p NUT LENGTHS ) -1 ) 

[33] PATH LTH-OUTPUT i NPOINTS+ 1 ; 7 ] 

[34] -0FFLINE*0PRINT>2 

[35] OFFLINE :CH-NUT2CH 

[36] 0 PRINT NUT2FMT FMT OUTPUT 

V 
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VDISTANCEl E]]V 
7 DIST+DI STANCE POINTS 

[1] fl THIS FUNCTION GENERATES AN N*N DISTANCE MATRIX 

[2] DISTH POINT Si ;2]° . -POINT Si ; 2 ] ) *2 

[ 3 ] DIST+-DIST+ ( POINT Si ; 1 ] * . -POINT Si ; 1 ] ) *2 

[4] DIST+DIST* 0.5 

V 
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